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Abstract 

This  report  reviews  the  initial  snowpack  distribution  assumed  to  be  in 
place  and  available  for  melt  in  the  probable  maximum  flood  (PMF)  analy¬ 
sis  of  Success  Dam  conducted  by  the  Corps  of  Engineers.  This  project  re¬ 
viewed  two  aspects  of  the  initial  snow  distribution:  the  snow-covered  area 
(SCA)  distribution  and  the  snow  water  equivalent  (SWE)  distribution.  We 
analyzed  satellite  imagery  to  determine  the  daily  SCA  of  the  Success  Dam 
watershed  from  2000  to  2014.  This  analysis  was  based  on  time-domain  fil¬ 
tering  of  NASA’s  daily  snow  product.  The  SCA  of  the  Success  Dam  water¬ 
shed  is  highly  dynamic  with  large  day-to-day  variations.  The  maximum 
95th  percentile  SCA  for  each  elevation  band  was  selected  to  describe  the 
initial  snow  cover  in  place  at  the  start  of  the  PMF  analysis.  We  next  re¬ 
viewed  the  initial  SWE  distribution.  Above  5000  ft,  the  initial  SWE  distri¬ 
bution  was  greater  than  the  snowmelt  that  occurred  during  the  PMF  and 
was  therefore  heat  transfer  limited.  In  this  region  the  snowmelt  runoff 
during  the  PMF  was  insensitive  to  the  amount  of  SWE.  Below  5000  ft,  the 
snowmelt  was  limited  by  the  initial  SWE.  Runoff  from  this  region  is  sensi¬ 
tive  to  the  amount  of  initial  SWE. 
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Executive  Summary 

This  report  reviews  the  initial  snowpack  distribution  assumed  to  be  in 
place  and  available  for  melt  in  the  probable  maximum  flood  (PMF)  analy¬ 
sis  of  Success  Dam  conducted  by  the  U.S.  Army  Corps  of  Engineers 
(USACE  2013).  We  conducted  our  review  because  of  the  significant  impact 
of  snowmelt  on  the  PMF  results  and  the  scarcity  of  data  available  to  de¬ 
scribe  the  snowpack  conditions.  We  reviewed  two  aspects  of  the  initial 
snow  distribution:  the  snow-covered  area  (SCA)  distribution  and  the 
snow-water  equivalent  (SWE)  distribution  over  the  Success  Dam  water¬ 
shed. 

The  Corps  PMF  snowmelt  analysis  (USACE  2013)  divided  the  watershed 
into  ten  1000  ft  elevation  bands.  The  PMF  assumed  that  an  initial  SWE 
depth  was  in  place  in  each  elevation  band  and  was  available  to  melt 
(USACE  2013).  The  initial  SWE  ranged  from  a  maximum  of  54.4  in.  at 
8000-9000  ft  to  a  minimum  of  4  in.  at  3000-4000  ft.  The  Corps  PMF 
snowmelt  analysis  assumed  that  no  snow  was  in  place  below  3000  ft.  Each 
elevation  band  was  implicitly  assumed  to  be  100%  snow  covered.  The  ini¬ 
tial  distribution  of  the  snowpack  resulted  in  232.7  square  miles  (59.7%)  of 
the  watershed  covered  with  snow  and  a  total  SWE  volume  of  267  x  ios 
acre-ft. 

For  our  review,  we  at  the  U.S.  Army  Cold  Regions  Research  and  Engineer¬ 
ing  Laboratory  (CRREL)  analyzed  satellite  imagery  to  determine  the  daily 
SCA  of  the  Success  Dam  watershed  from  2000  to  2014.  This  analysis  was 
based  on  time-domain  filtering  (TDF)  of  NASA’s  (National  Aeronautics 
and  Space  Administration)  daily  MOD10A1  snow  product.  TDF  results  in 
virtually  cloud-free  daily  snow  maps  displaying  the  fractional  SCA  (fSCA) 
of  pixels  not  obscured  by  clouds.  The  pixel  size  is  roughly  500  x  500  m, 
and  the  Success  Dam  watershed  is  covered  by  4710  pixels.  We  estimated 
the  fSCA  of  each  pixel  for  each  day  for  2000-2014. 

The  CRREL  SCA  analysis  found  that  the  SCA  of  the  Success  Dam  water¬ 
shed  is  highly  dynamic  with  large  day-to-day  variations.  This  study  consid¬ 
ered  various  options  for  selecting  the  appropriate  SCA  to  use  in  determin¬ 
ing  the  initial  snow  distribution.  We  calculated  the  SCA  statistics  for  each 
day  of  the  winter  season  and  selected  the  maximum  95th  percentile  SCA 
for  each  elevation  band  to  describe  the  initial  snow  cover  in  place  at  the 


ERDC/CRREL  TR-15-16 


ix 


start  of  the  analysis  as  being  sufficiently  conservative  without  being  ex¬ 
treme.  The  total  watershed  SCA  that  resulted  from  selecting  the  maximum 
95th  percentile  SCA  for  each  elevation  band  was  101.6  square  miles 
(26.03%  °f  the  watershed).  This  is  greater  than  the  total  watershed  SCA 
observed  on  99.5%  of  the  days  with  SCA  greater  than  zero.  The  SWE  vol¬ 
ume  that  resulted  from  combining  the  SWE  depth  for  each  elevation  band 
used  in  the  PMF  study  with  the  95th  percentile  SCA  for  each  elevation 
band  was  178  x  io3  acre-ft. 

We  next  reviewed  the  initial  SWE  distribution  used  in  the  Corps  PMF 
study  (USACE  2013).  Above  5000  ft,  the  initial  SWE  distribution  was 
greater  than  the  snowmelt  that  occurred  during  the  PMF.  The  snowmelt  in 
this  region  was  therefore  heat  transfer  limited.  In  this  region,  the  snow¬ 
melt  runoff  during  the  PMF  was  insensitive  to  the  amount  of  SWE  given 
that  the  SWE  was  greater  than  the  snowmelt  amount.  Below  5000  ft,  the 
snowmelt  was  limited  by  the  initial  SWE.  Runoff  from  this  region  was  sen¬ 
sitive  to  the  amount  of  initial  SWE.  Unfortunately  there  is  no  information 
on  the  actual  SWE  in  this  region  of  the  Success  Dam  watershed.  Based  on 
the  work  of  others  and  the  SCA  analysis  described  in  this  report,  the  distri¬ 
bution  selected  in  this  region  seems  appropriate. 

The  100%  snow  cover  area  used  in  the  PMF  study  is  also  a  requirement  of 
the  HEC-HMS  elevation-band  snow  model  (USACE  2010a).  HEC-HMS  is 
not  currently  capable  of  including  the  impact  of  SCA  in  the  elevation-band 
snow-model  option  that  was  used  for  the  PMF  runoff  calculations.  An  al¬ 
ternative  is  to  use  the  gridded  snow-model  option  instead. 
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1  Introduction 

Success  Dam,  a  U.S.  Army  Corps  of  Engineers  (USACE)  flood  control  pro¬ 
ject,  is  located  in  southern  California  on  the  Tule  River.  This  report  re¬ 
views  the  initial  snowpack  distribution  assumed  to  be  in  place  and  availa¬ 
ble  for  melt  in  the  probable  maximum  flood  (PMF)  analysis  of  Success 
Dam  conducted  by  the  Corps  of  Engineers  (USACE  2013). 

A  significant  amount  of  the  PMF  runoff  from  the  watershed  upstream  of 
Success  Dam  is  from  snowmelt  (USACE  2013).  Therefore,  we  conducted 
our  review  because  of  the  significant  impact  of  snowmelt  on  the  PMF  re¬ 
sults  and  the  scarcity  of  data  available  to  describe  the  snowpack  condi¬ 
tions.  The  Corps  PMF  snowmelt  analysis  divided  the  watershed  into  ten 
1000  ft  elevation  bands,  assigning  an  initial  snow-water  equivalent  (SWE) 
depth  to  each  band.  The  initial  SWE  ranged  from  a  maximum  54.4  in.  at 
8000-9000  feet  to  a  minimum  of  4  in.  at  3000-4000  ft.  The  Corps  PMF 
snowmelt  analysis  assumed  that  no  snow  was  in  place  below  3000  ft.  Each 
elevation  band  was  implicitly  assumed  to  be  100%  snow  covered.  The  ini¬ 
tial  distribution  of  the  snowpack  resulted  in  232.7  square  miles  (59.7%)  of 
the  watershed  covered  with  snow  and  a  total  SWE  volume  of  267  x  ios 
acre-ft. 

The  initial  snowpack  distribution  resulted  from  the  combination  of  SWE 
and  snow-cover  area  (SCA)  assigned  to  each  elevation  band.  The  SWE  as¬ 
signed  to  each  elevation  band  in  the  PMF  analysis  of  Success  Dam  does 
not  seem  unreasonable.  However,  the  implicit  assumption  of  100%  SCA  in 
each  elevation  band  is  likely  to  be  unduly  conservative.  The  impact  of  SCA 
on  snowmelt  calculations  has  long  been  recognized  (see,  for  example,  An¬ 
derson  1973).  However,  it  has  not  been  practical  to  quantitatively  observe 
the  SCA  of  a  moderate-to-large-size  watershed  at  a  resolution  and  fre¬ 
quency  that  is  useful.  As  a  result,  the  snowmelt  component  of  the  HEC- 
HMS  elevation-band  snow  model  assumes  that  each  subbasin  or  elevation 
band  has  100%  SCA  at  all  times  (USACE  2010a). 

A  new  satellite  imagery  analysis  tool  developed  at  the  U.S.  Army  Cold  re¬ 
gions  Research  and  Engineering  Laboratory  (CRREL)  for  application  in  re¬ 
mote  austere  watersheds  now  allows  the  SCA  to  be  automatically  mapped 
using  NASA’s  (National  Aeronautics  and  Space  Administration)  daily 
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MOD10A1  snow  product.  This  process,  referred  to  as  Time-Domain  Filter¬ 
ing  (TDF),  results  in  virtually  cloud-free  daily  snow  maps.  The  daily  snow 
maps  display  the  fractional  SCA  (fSCA)  at  a  pixel  resolution  of  500  x  500 
m,  and  the  Success  Dam  watershed  is  covered  by  4710  pixels.  We  used 
TDF  to  determine  the  fSCA  of  each  pixel  for  each  day  for  2000-2014. 

The  CRREL  SCA  analysis  found  that  the  SCA  of  the  Success  Dam  water¬ 
shed  is  highly  dynamic  and  varies  significantly  from  day  to  day  but  that 
the  SCA  analysis  was  consistent  with  the  topography  of  the  watershed  and 
the  estimated  air  temperature  range.  We  considered  various  options  for 
selecting  the  appropriate  SCA  to  use  in  determining  the  initial  snowpack 
distribution.  The  SCA  selected  is  conservative  but  not  extreme,  and  it  sig¬ 
nificantly  modifies  the  snowpack  distribution  compared  to  the  initial 
snowpack  distribution  assumed  to  be  in  place  at  the  start  of  the  PMF  anal¬ 
ysis  of  Success  Dam  (USACE  2013). 

We  next  reviewed  the  initial  SWE  distribution  used  in  the  Corps  PMF 
study  (USACE  2013).  Above  5000  ft,  the  initial  SWE  distribution  was 
greater  than  the  snowmelt  that  occurred  during  the  PMF.  The  snowmelt  in 
this  region  was  therefore  heat  transfer  limited.  In  this  region,  the  snow¬ 
melt  runoff  during  the  PMF  was  insensitive  to  the  amount  of  SWE  given 
that  the  SWE  was  greater  than  the  snowmelt  amount.  Below  5000  ft,  the 
snowmelt  was  limited  by  the  initial  SWE.  Runoff  from  this  region  was  sen¬ 
sitive  to  the  amount  of  initial  SWE.  Unfortunately  there  is  no  information 
on  the  actual  SWE  in  this  region  of  the  Success  Dam  watershed.  Based  on 
the  work  of  others  and  the  SCA  analysis  described  in  this  report,  the  distri¬ 
bution  selected  in  this  region  seems  appropriate. 
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2  Success  Dam  Watershed 

The  Success  Dam  watershed  is  located  on  the  western  slope  of  the  south¬ 
ern  Sierra  Nevada  Mountains  in  California  (USACE  2013).  The  390- 
square-mile  drainage  area  above  the  dam  is  largely  steep  mountainous  ter¬ 
rain.  Watershed  elevations  range  from  about  10,000  ft  in  the  mountains  to 
550  ft  at  the  dam.  Figure  1  shows  the  layout  of  the  watershed;  and  the 
USACE  reports  (1998,  2013)  further  describe  the  watershed  characteris¬ 
tics,  vegetation,  climate,  precipitation,  and  flooding.  The  locations  of  data 
collection  stations  in  the  watershed  are  shown  in  Figure  1  and  are  de¬ 
scribed  in  Table  1. 


Figure  1.  The  Lake  Success  watershed. 
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Table  1.  SWE  and  Air  Temperature  stations  located  in  Success  Dam  watershed. 


ID 

Station 

Parameter 

Rate 

Longitude 

Latitude 

Elevation 

(ft) 

Period  of 
Record 

Operator 

OEM 

Old  Enterprise 
Mill 

SWE 

Yearly/ 

Monthly 

-118.678 

36.243 

6600 

1937-2014 

CA  Dept,  of 
Forestry  and  Fire 
Protection 

QKA 

Quaking  Aspen 

SWE 

Yearly/ 

Monthly 

-118.545 

36.122 

7000 

1937-2014 

Tule  River 

Ranger  District 

QUA 

Quaking  Aspen 

SWE 

Daily 

-118.53994 

36.11754 

7200 

1981-2014 

CA  Dept,  of 

Water  Resources 

EGL 

Eagle  Creek 

Air  Temp. 

Hourly 

-118.641 

35.983 

6700 

1995-2009 

USACE 

HSS 

Hossack 

Air  Temp. 

Hourly 

-118.619 

36.181 

7100 

1995-2009 

USACE 

MNH 

Mountain 

Home 

Air  Temp. 

Hourly 

-118.71 

36.242 

5400 

1995-2009 

USACE 

RGC 

Rogers  Camp 

Air  Temp. 

Hourly 

-118.637 

36.108 

6200 

1995-2009 

USACE 

sew 

Success 

AirTemp. 

Hourly 

-118.9139 

36.05 

590 

1995-2009 

USACE 

2.1  Topography 

We  analyzed  the  topography  of  the  Success  Dam  watershed  on  a  pixel-by¬ 
pixel  basis  using  pixels  defined  in  the  MOD10A1  snow  product  imagery. 
MODioAl  snow  product  imagery  is  in  the  form  of  “tiles”  with  each  tile  ap¬ 
proximately  1200  x  1200  km  in  area  and  containing  approximately  2400  x 
2400  pixels.  Each  pixel  is  roughly  500  x  500  m  in  size.  The  Sacramento 
District  provided  a  shapefile  of  the  Success  Dam  watershed  boundary. 
Combining  the  watershed  shapefile  and  the  MODioAl  tile  produced  a 
shapefile  that  contained  4710  polygons  that  described  each  MODIS  pixel 
within  the  watershed.  Only  the  MODIS  pixels  whose  center  point  fell 
within  the  watershed  were  included  in  the  analysis.  We  used  ArcMap  10. 1 
SPi  to  determine  the  elevation,  slope,  and  aspect  of  each  MODIS  pixel  by 
using  a  DEM  called  "mos_prj_ft"  that  the  Sacramento  District  provided 
(A.  Duren,  personal  communication).  The  cell  size  of  this  DEM  was  a  little 
over  30  ft.  We  determined  the  elevation,  slope,  and  aspect  of  each  pixel  us¬ 
ing  ArcGIS  (ESRI  2011).  Specifically,  we  calculated  the  elevation  and  slope 
by  using  the  “Zonal  Statistics  as  Table”  tool;  and  we  calculated  the  mean 
aspect  by  using  a  custom  script  tool  (Beyerhelm  2013). 

Figure  2  shows  the  distribution  of  the  pixel  aspects.  The  watershed  faces 
largely  to  the  southwest.  We  also  followed  the  example  of  the  PMF  study 
(USACE  2013)  and  divided  the  watershed  into  1000  ft  elevation  bands. 
Figure  3  and  Table  2  show  the  area  of  each  elevation  band. 
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Figure  2.  The  aspect  distribution  of  Success  Dam  watershed  pixels. 
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Figure  3.  The  area  of  each  elevation  band. 
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Table  2.  The  area  of  each  elevation  band  and  the  total  area. 


Elevation  Band 

Area 

(Square 

Miles) 

Cumulative  Area 
From  the  Peak 
(Square  Miles) 

515-1000 

31.9 

390.1 

1000-2000 

66.7 

358.2 

2000-3000 

58.8 

291.5 

3000-4000 

44.2 

232.7 

4000-5000 

52.0 

188.5 

5000-6000 

51.9 

136.5 

6000-7000 

40.3 

84.6 

7000-8000 

26.0 

44.3 

8000-9000 

14.2 

18.2 

9000-10,000 

4.1 

4.1 

2.2  Air  temperature 

We  calculated  the  statistics  of  the  daily  average  air  temperature  for  each 
day  of  the  winter  season  to  provide  an  overview  of  the  wintertime  condi¬ 
tions  by  using  the  “Cyclic  Analysis  Statistics  Math  Tool”  of  HEC-DSSVue 
(USACE  2009).  Hourly  air  temperature  records  are  available  at  five  sta¬ 
tions  in  the  Success  Dam  watershed  (Table  l ).  The  statistics  were  calcu¬ 
lated  based  on  the  entire  period  of  record  for  each  station  listed  in  Table  l. 
The  daily  average  air  temperature  was  rarely  lower  than  32°F  (indicated  in 
Figure  4)  at  any  of  the  stations  throughout  the  winter.  It  is  interesting  to 
note  the  relatively  rapid  rise  of  the  temperature  at  the  Success  Dam  com¬ 
pared  to  the  other  stations  starting  in  mid-winter  and  continuing  through 
early  summer.  Figure  5  shows  the  daily  average  minimum  temperature  ob¬ 
served  at  each  station  for  each  day  of  the  winter  season.  The  daily  average 
minimum  air  temperatures  were  below  32°F  for  all  the  stations  except 
Success  Dam  for  December  through  April.  These  higher  temperatures  at 
Success  Dam  attest  to  the  relatively  warm  conditions  at  the  downstream 
and  lowest  end  of  the  watershed. 

The  hourly  elevation  of  the  32°F  isotherm  of  the  watershed  was  estimated 
through  linear  interpolation  of  the  hourly  air  temperatures  recorded  at  the 
five  stations.  The  elevation  of  the  stations  ranged  from  590  to  7200  ft  (alt¬ 
hough  the  highest  four  stations  are  all  located  above  5400  ft).  We  found 
that  elevations  above  the  32°F  isotherm  were  always  colder  than  eleva¬ 
tions  below  and  that  temperature  inversions  were  all  short  lived;  none 
lasted  longer  than  a  day.  The  32°F  isotherm  provides  information  on  the 
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conditions  that  control  heat  transfer  between  the  snowpack  and  the  at¬ 
mosphere  and  the  snowpack  persistence  as  snow  is  more  likely  to  melt 
when  exposed  to  air  temperatures  above  32°F.  Figure  6  shows  the  mini¬ 
mum,  maximum,  and  average  elevations  of  the  daily  32°F  isotherm 
throughout  the  winter.  The  average  32°F  isotherm  rarely  falls  below  6000 
ft,  while  the  minimum  has  reached  as  low  as  2000  ft  in  isolated  cases. 


Figure  4.  Daily  average  air  temperatures  for  each  day  of  the  winter. 


Figure  5.  Minimum  daily  average  air  temperatures  for  each  day  of  the  winter. 
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Figure  6.  The  elevation  of  the  32°F  isotherm. 
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3  PMF  Initial  Snow  Conditions 


The  initial  snowpack  distribution  assumed  to  be  in  place  at  the  start  of  the 
PMF  analysis  of  Success  Dam  (USACE  2013)  are  shown  in  Table  3  and  dis¬ 
played  in  Figure  7.  In  that  initial  analysis,  the  volume  of  SWE  for  each  ele¬ 
vation  band  was  found  by  multiplying  the  SWE  assigned  to  each  band 
(Table  3)  by  the  area  of  the  elevation  band  (Table  2). 


Table  3.  PMF  Starting  Conditions  (USACE  2013). 


Project  Datum 
(ft) 

Snow  Density 
(ft) 

Snow  Depth 
(in.) 

Snow-Water  Equivalent 
(in.) 

515-1000 

No  Snow  on  the  Ground 

0 

1000-2000 

No  Snow  on  the  Ground 

0 

2000-3000 

No  Snow  on  the  Ground 

0 

3000-4000 

0.40 

10 

4.0 

4000-5000 

0.40 

20 

8.0 

5000-6000 

0.39 

40 

15.6 

6000-7000 

0.37 

90 

33.3 

7000-8000 

0.35 

140 

49.0 

8000-9000 

0.32 

170 

54.4 

9000-10,000 

0.30 

175 

52.5 

Figure  7.  The  PMF  SWE  distribution  by  elevation. 
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4  Snow-Covered-Area  Analysis 

4.1  Background 

Our  analysis  of  SCA  is  based  on  the  NASA-developed  SCA  product, 
MODioAl,  the  MODIS/Terra  Snow  Cover  Daily  L3  Global  500  m  Grid 
(Hall  et  al.  2006),  which  provides  a  daily  estimate  of  the  fSCA  of  the  globe 
at  a  500  m  pixel  resolution.  The  SCA  is  estimated  using  the  normalized 
difference  snow  index  (NDSI),  which  leverages  the  contrasting  reflectance 
of  snow  in  the  visible  and  shortwave-infrared  regions  of  the  spectrum 
(Hall  et  al.  2006).  The  fSCA  is  estimated  using  an  empirical  relationship 
between  fSCA  and  NDSI.  This  process  may  have  difficulties  in  heavily  for¬ 
ested  areas  (Salomonson  and  Appel  2004).  Areas  obscured  by  clouds  are 
masked  out  of  the  MODioAl  product  based  on  the  MODIS  cloud  mask 
(MOD35).  The  MODIS  cloud  mask  (MOD35)  uses  a  host  of  algorithms  and 
19  of  MODIS’s  36  bands  to  determine  if  clouds  are  contained  in  the 
MODIS  image  (Ackerman  et  al.  2010). 

The  large  and  changeable  areas  obscured  by  clouds  in  each  daily  SCA 
product  represent  a  significant  impediment  to  the  use  of  the  MODioAl 
SCA  product  in  water  resources  investigations.  CRREL  has  developed  a 
new  technique,  referred  to  as  Time  Domain  Filtering  (TDF),  to  address 
this  issue  of  cloud  obscuration  (Morriss  et  al.,  forthcoming).  In  general,  we 
find  mean  cloud  persistence,  or  the  average  time  between  cloud-free 
scenes,  to  be  less  than  a  day  at  Terra  MODIS’s  temporal  scale  (approxi¬ 
mately  1  day).  This  cloud-cover  transience  gives  us  the  opportunity  to  fill 
in  gaps  in  the  current  day’s  imagery  with  recent  data.  Because  we  use  mul¬ 
tiple  days’  imagery  as  inputs,  we  choose  to  limit  our  outputs  to  SCA  that  is 
persistent  on  a  multi-day  scale.  Each  output  pixel’s  value  is  based  on  up  to 
16  daily  images  though  only  the  two  most  recent  cloud-free  images  influ¬ 
ence  the  value  of  the  outputs.  This  approach  limits  errors  of  commission  in 
snow  cover  that  arise  from  cloud  and  cloud  shadow  misclassification  in  the 
MODioAl  detection  algorithms.  If  SCA  is  present  in  both  of  the  two  most 
recent  cloud-free  images,  the  output  pixel’s  value  is  set  to  the  fSCA  value  of 
the  most  recent  image.  If  either  pixel’s  value  in  the  cloud-free  pair  is  snow 
free,  the  output  pixel  is  assumed  to  be  snow  free.  This  requirement  some¬ 
times  filters  transient  SCA  from  the  output  product  in  addition  to  spurious 
SCA  but  creates  a  conservative  and  robust  delineation  of  SCA  extent.  In 
the  rare  occurrence  that  a  pixel  does  not  have  any  cloud-free  readings  in 
the  16  input  scenes,  the  output  is  cloud. 
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4.2  Overview 

Our  study  analyzed  SCA  within  the  Success  Dam  watershed  by  determin¬ 
ing  the  fSCA  for  each  pixel  for  each  day  by  using  the  TDF  analysis  of  the 
MOD10A  snow  product.  We  determined  the  fSCA  for  each  of  the  4710 
MODIS  pixels  in  the  watershed  for  each  of  the  5405  days  in  the  period  of 
record,  which  spanned  from  4  March  2000  through  20  December  2014. 
There  were  45  days  of  fSCA  data  missing  due  to  gaps  in  the  MOD10A  data. 
There  were  2715  days  with  SCA  of  which  2133  days  (78.6%)  had  unique 
data;  582  days  with  SCA  (21.4%)  were  duplicates  caused  by  the  presence  of 
clouds. 

We  determined  the  total  area  covered  by  SCA  for  each  day  of  the  period  of 
record  by  summing  the  fSCA  of  all  the  pixels  in  the  watershed  and  multi¬ 
plying  the  result  by  the  area  of  a  MODIS  pixel.  Figure  8  shows  the  time  se¬ 
ries  of  total  watershed  SCA.  It  is  interesting  to  note  the  “spikiness”  of  the 
total  SCA.  The  area  of  the  watershed  that  is  continuously  snow  covered 
each  winter  is  relatively  small,  generally  less  than  20  square  miles.  Periods 
when  much  larger  areas  are  snow  covered  last  only  a  few  days. 


Figure  8.  The  total  watershed  SCA. 


Figure  9  shows  the  time  series  of  total  watershed  SCA  as  a  percentage  of 
the  total  Success  Dam  watershed.  This  demonstrates  that  the  SCA  gener¬ 
ally  ranges  between  5%  and  10%  of  the  watershed  area.  Figure  10  displays 
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the  statistics  of  SCA  for  each  day  in  the  winter  and  includes  the  daily  aver¬ 
age  SCA  and  the  maximum  and  minimum  recorded  SCA. 


Figure  9.  The  total  watershed  SCA  as  a  percentage  of  the  watershed  area. 


Figure  10.  Average,  maximum,  and  minimum  total  watershed  SCA  for  each  day  of 

the  winter. 
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4.3  Daily  SCA  results 

For  demonstrational  purposes,  we  selected  several  days  in  the  period  of 
record  to  view  in  detail.  On  29  January  2001,  the  largest  total  watershed 
SCA  over  the  period  of  record  occurred  (Figure  11).  A  second  image  dis¬ 
plays  the  SCA  map  seven  days  later  on  5  February  2001  (Figure  11).  The 
dynamic  nature  of  the  SCA  can  be  seen  in  the  dramatic  reduction  in  SCA 
over  these  seven  days. 


Figure  11.  The  SCA  on  29  January  2001  and  seven  days  later  on  5  February  2001. 


Figure  12  shows  the  estimated  SCA  for  12  January  2005.  This  is  also  a  date 
when  one  of  the  largest  watershed  SCA  was  recorded.  The  pixels  marked 
as  cloud  covered  mean  that  the  MODIS  cloud  mask  interpreted  the  condi¬ 
tions  over  these  pixels  as  cloud  covered  for  16  consecutive  days  or  more  on 
12  February  2005.  Note  that  some  areas  at  higher  elevations  are  snow  free 
and  that  these  areas  match  snow  free  areas  recorded  on  29  January  2001 
(Figure  11). 

Figure  13  shows  the  estimated  SCA  for  18  February.  This  date  also  had  one 
of  the  largest  watershed  SCA  recorded.  In  large  part,  the  SCA  follows  the 
elevation  contours,  but  not  exclusively  so.  There  are  several  instances 
where  the  SCA  extends  downhill  in  isolated  salients,  essentially  cutting 
across  the  elevation  contours.  These  downhill  SCA  projections  may  result 
from  spatial  variations  in  the  distribution  of  snowfall  that  are  not  entirely 
coincidental  with  the  contour  locations. 
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Figure  12.  The  SCA  on  12  January  2005. 
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Figure  13.  The  SCA  on  18  February  2009. 
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4.4  General  SCA  distribution 

We  first  analyzed  the  general  distribution  of  SCA  throughout  the  water¬ 
shed  by  calculating  the  longest  duration  each  winter  that  each  pixel  was 
continuously  covered  by  snow  (fSCA  >  20%).  We  then  averaged  these  du¬ 
rations  over  the  period  of  record  and  mapped  them  (Figure  14).  This  map 
also  includes  the  elevation  contours,  showing  that  the  distribution  of  SCA 
very  closely  follows  these  contours.  This  close  agreement  of  the  SCA  dura¬ 
tion  and  the  contours  is  taken  as  a  partial  confirmation  of  the  TDF  pro¬ 
cess. 


We  used  the  total  watershed  average  daily  SCA  shown  in  Figure  10  to  esti¬ 
mate  the  elevation  of  the  average  effective  daily  snow  line.  Here,  the  effec¬ 
tive  snow  line  is  defined  by  the  elevation  above  which  the  area  equals  the 
total  watershed  SCA.  We  used  the  information  in  column  3  of  Table  2  to 
interpolate  the  elevation  that  corresponded  to  the  area  of  the  total  water¬ 
shed  SCA  of  each  day.  This  is  a  rough  measure  of  the  snow  line,  but  the 
close  relation  of  the  snow  line  and  the  32°F  isotherm  can  be  seen  in  Figure 
15.  These  close  results  also  suggest  that  the  TDF  process  provided  reasona¬ 
ble  and  accurate  results. 


Figure  14.  The  average  continuous  SCA  duration  each  winter  by  pixel. 
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Figure  15.  The  average  elevation  of  the  32°F  isotherm  and  the  snow  line. 
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We  also  calculated  the  SCA  for  each  elevation  band  for  each  day  over  the 
period  of  record,  having  determined  the  elevation  for  each  pixel  as  de¬ 
scribed  earlier.  The  SCA  for  an  elevation  band  on  day  i,  SCA'eb ,  was 


i 


(1) 


where 


n  =  the  number  of  pixels  in  the  elevation  band  and 
fSCAj  =  the  fSCA  of  the  /th  pixel. 

We  then  determined  the  SCA  statistics  for  each  day  of  the  year  by  using 
the  “Cyclic  Analysis  Statistics  Math  Tool”  of  HEC-DSSVue  (USACE  2009). 
Figure  16  shows  the  daily  average  SCA  for  each  elevation  band.  In  this  fig¬ 
ure,  the  daily  average  SCA  for  each  elevation  band  drops  rapidly  as  the  ele¬ 
vation  band  elevation  decreases,  in  agreement  with  Figure  14.  Figure  17 
shows  the  daily  95th  percentile  SCA  for  each  elevation  band.  The  95th  per¬ 
centile  SCA  is  greater  than  95%  of  the  SCA  recorded  on  that  day  of  the 
year.  The  95th  percentile  SCA,  while  greater  in  value  than  the  average,  also 
drops  rapidly  as  the  band  elevation  decreases,  especially  below  7000  ft. 
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Figure  16.  The  average  SCA  on  each  day  of  winter  for  each  elevation  band. 
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Figure  17.  The  95th  percentile  SCA  on  each  day  of  the  winter  for  each  elevation  band. 
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5  Application  of  SCA  to  the  PMF 
Determination 

The  starting  conditions  of  the  PMF  study  of  the  Success  Dam  (USACE 
2013)  had  snow  extending  from  the  highest  elevations  down  to  3000  ft  as 
shown  in  Table  3.  There  was  no  accounting  for  the  likely  snow  covered 
area  in  each  band.  The  SCA  for  each  band  was  implicitly  assumed  to  be 
100%.  As  can  be  seen  in  the  previous  section,  the  likelihood  of  SCA  de¬ 
creases  rapidly  as  the  elevation  decreases. 

There  is  no  discussion  regarding  the  inclusion  of  SCA  in  a  PMF  study  in 
any  of  the  relevant  guidance  documents  of  the  USACE.  Therefore,  we  de¬ 
cided  to  review  four  possible  options  for  selecting  the  SCA  for  the  elevation 
bands.  The  results  are  shown  in  Table  4  and  are  displayed  in  Figure  18  and 
Figure  19. 

5.1  SCENARIO  A 

We  based  this  scenario  on  the  average  SCA  in  each  band  that  was  observed 
on  the  day  of  the  maximum  SWE  observed  at  the  Quaking  Aspen  station 
(QUA)  over  the  period  from  2000  to  2014.  As  it  is  likely  that  extreme 
events  will  happen  at  or  around  the  time  of  maximum  SWE,  we  selected 
the  SCA  corresponding  to  the  date  of  the  maximum  SWE  for  each  eleva¬ 
tion  band  for  each  year  from  2000  through  2014.  We  averaged  the  annual 
SCAs  to  arrive  at  a  final  value  for  each  elevation  band.  These  results  were 
not  particularly  conservative  as  the  total  watershed  SCA  of  this  scenario  is 
only  greater  than  the  total  watershed  SCA  observed  on  82%  of  the  days 
with  a  total  SCA  greater  than  zero. 

5.2  SCENARIO  B 

We  based  this  scenario  on  the  maximum  SCA  in  each  band  that  was  ob¬ 
served  on  the  day  of  the  maximum  SWE  observed  at  the  Quaking  Aspen 
station  (QUA)  for  the  years  2000  through  2014.  This  is  similar  to 
SCENARIO  A  except  that  we  selected  the  maximum  SCA  observed  on  the 
date  of  the  maximum  SWE.  The  results  are  shown  in  Table  4  and  are  dis¬ 
played  in  Figure  18.  These  results  are  more  conservative  than  SCENARIO 
A.  The  total  watershed  SCA  of  this  scenario  is  greater  than  the  total  water¬ 
shed  SCA  observed  on  96.3%  of  the  days  with  a  total  SCA  greater  than 


zero. 
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5.3  SCENARIO  C 

This  scenario  used  the  average  maximum  30-day  SCA.  The  daily  SCA  for 
each  elevation  band  was  converted  to  a  time  step  of  30  days  with  the  value 
of  each  30-day  variable  set  equal  to  the  average  of  the  30  days  that  it 
spanned.  Then  the  maximum  30-day  value  was  selected  for  each  winter. 
The  average  over  all  winters  was  then  calculated  to  arrive  at  the  average 
maximum  30-day  SCA.  These  results  are  more  conservative  than 
SCENARIO  A  but  slightly  less  than  SCENARIO  B.  The  total  watershed  SCA 
of  this  scenario  is  greater  than  the  total  watershed  SCA  observed  on  94.1% 
of  the  days  with  a  total  SCA  greater  than  zero. 

5.4  SCENARIO  D 

This  scenario  used  the  maximum  95th  percentile  SCA.  We  selected  the 
maximum  95th  percentile  SCA,  displayed  in  Figure  17,  from  each  elevation 
band.  These  results  are  more  conservative  than  the  other  scenarios.  The 
total  watershed  SCA  of  this  scenario  is  greater  than  the  total  watershed 
SCA  observed  on  99.5%  of  the  days  with  a  total  SCA  greater  than  zero. 

5.5  Selection  of  the  SCA  Scenario 

The  selection  of  the  SCA  scenario  to  use  is  complicated  because  of  the 
highly  dynamic  nature  of  SCA  in  the  Success  Dam  watershed  as  shown  in 
Figure  8  and  Figure  9.  The  accumulation  of  SWE  requires  time  for  succes¬ 
sive  snowfall  events  to  accumulate.  Highly  ephemeral  SCA  responds  to 
each  snowfall  event  but  quickly  ablates  before  significant  SWE  can  accu¬ 
mulate.  The  SCA  scenario  selected  has  to  strike  a  balance  between  being 
sufficiently  conservative  in  representing  the  largest  area  of  the  watershed 
that  could  be  covered  with  snow  when  the  PMF  occurred  and  having  suffi¬ 
cient  time  in  place  for  the  SWE  to  accumulate.  We  selected  the  maximum 
95th  percentile  SCA  (SCENARIO  D)  for  each  elevation  band  as  being  suffi¬ 
ciently  conservative  without  representing  transient  events  where  some  of 
the  SCA  had  little  SWE.  The  total  watershed  SCA  that  results  from  select¬ 
ing  the  maximum  95th  percentile  SCA  for  each  elevation  band  is  101.6 
square  miles  (26.03%  of  the  watershed).  This  is  greater  than  the  total  wa¬ 
tershed  SCA  observed  on  99.5%  of  the  days  with  SCA  greater  than  zero. 

The  SWE  volume  that  results  from  combining  the  SWE  depth  used  in  the 
PMF  study  with  the  95th  percentile  SCA  for  each  elevation  band  is  178  x 
103  acre-ft. 
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Table  4.  The  SCA  by  elevation  and  results  for  each  scenario. 


Elevation  Band 

PMF 

(%) 

A 

(%) 

B 

(%) 

c 

(%) 

D 

(%) 

515-1000 

- 

0.0 

0 

0.0 

0.0 

1000-2000 

- 

0.0 

0 

0.0 

0.0 

2000-3000 

- 

0.0 

0 

0.0 

0.8 

3000-4000 

100 

0.0 

0.5 

2.9 

7.7 

4000-5000 

100 

0.1 

0.5 

6.8 

22.0 

5000-6000 

100 

1.6 

9.7 

11.2 

40.7 

6000-7000 

100 

18.0 

45.4 

27.4 

60.6 

7000-8000 

100 

54.0 

83.0 

64.9 

87.6 

8000-9000 

100 

74.8 

97.4 

84.8 

97.7 

9000-10,000 

100 

87.7 

99.3 

93.3 

99.0 

PMF 

A 

B 

c 

D 

Total  SCA 
(square  miles) 

232.7 

36.4 

63.3 

54.4 

101.6 

Total  SCA 

(%) 

59.7 

9.3 

16.2 

13.9 

26.0 

%  Days  Less 

Than  or  Equal  to 
Total  SCA 

100 

82.0 

96.3 

94.1 

99.5 

Total  SWE 
(acre-ft  x  io3) 

266.94 

91.09 

144.69 

115.87 

177.63 

Figure  18.  The  SCA  by  elevation  band  for  each  scenario. 


Elevation  Band 
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Figure  19.  The  volume  of  SWE  by  elevation  band  for  each  scenario. 
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6  SWE  Analysis 

This  section  reviews  the  initial  SWE  distribution  used  in  the  PMF  analysis. 
First,  we  compare  the  initial  SWE  to  the  SWE  measured  in  the  Success 
Dam  watershed  at  the  stations  listed  in  Table  l.  We  then  determine  the  im¬ 
pact  of  the  initial  SWE  distribution  on  the  snowmelt  calculated  during  the 
PMF. 

There  are  three  locations  where  SWE  is  observed  in  the  Success  Dam  wa¬ 
tershed  (Table  1):  an  automated  SWE  gage  providing  daily  observations 
(Quaking  Aspen  [QUA])  that  has  been  in  operation  since  1981  and  two 
snow  courses  (Quaking  Aspen  [QKA]  and  Old  Enterprise  Mill  [OEM]) 
where  data  have  been  collected  since  1937.  The  snow  course  data  were  col¬ 
lected  only  a  few  times  over  the  course  of  each  winter,  generally  around 
the  first  of  the  month. 

The  annual  maximum  SWE  measured  at  the  SWE  gage  (QUA)  and  the  two 
snow  courses  (QKA  and  OEM)  was  compiled  over  the  period  of  record  for 
each  source.  Observations  at  the  snow  courses  were  not  as  frequent  as  at 
the  snow  gage.  As  a  result,  the  snow  courses  did  not  capture  the  annual 
peak  as  faithfully  as  the  snow  gage.  The  average  annual  maximum  at  QKA 
is  about  93%  of  the  average  annual  maximum  recorded  at  QUA,  and  the 
average  annual  maximum  at  OEM  is  about  83%  of  the  average  maximum 
recorded  at  QUA.  We  should  note  that  the  all-time  maximum  SWE  rec¬ 
orded  by  any  of  the  three  sites  was  62.5  in.,  recorded  at  QKA  on  1  February 
1998. 

The  statistics  of  the  annual  maximum  SWE  were  calculated  using  a  log¬ 
normal  distribution  using  HEC-SSP  (USACE  2010b).  The  statistics  are 
listed  in  Table  5  and  shown  in  Figures  20-22.  Figure  23  shows  a  compari¬ 
son  of  the  SWE  statistics  and  the  initial  distribution  of  SWE  used  in  the 
PMF.  The  comparison  is  valid  only  for  the  elevation  bands  in  which  the 
SWE  sites  are  located.  It  can  be  seen  that  the  initial  SWE  falls  between 
about  a  10%  and  5%  annual  probability  (10-year  to  20-year  annual  return 
period).  This  comparison  in  and  of  itself  does  not  say  much  about  the  ap¬ 
propriateness  of  the  initial  SWE  distribution. 
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Table  5.  Statistics  of  annual  maximum  SWE  for  sites  in  the  Success  Dam  watershed. 


Annual  Probability 

Return  Period 

QKA 

OEM 

QUA 

(%) 

(years) 

(in.) 

(in.) 

(in.) 

0.2 

500 

83.1 

105.3 

109.3 

0.5 

200 

68.9 

85.9 

91.0 

1 

100 

59.0 

72.6 

78.3 

2 

50 

49.8 

60.4 

66.4 

5 

20 

38.6 

45.8 

51.8 

10 

10 

30.8 

35.8 

41.6 

20 

5 

23.5 

26.6 

31.8 

50 

2 

13.9 

15.1 

19.1 

80 

1.25 

8.2 

8.5 

11.5 

90 

1.11 

6.3 

6.4 

8.8 

95 

1.05 

5.0 

5.0 

7.1 

99 

1.01 

3.3 

3.1 

4.7 

Figure  20.  A  frequency  plot  for  Quaking  Aspen  (QUA). 

General  Frequency  Analytical  Plot  for  Success  Dam  Quaking  Aspen  CDEC 
Return  Period 


- 95  Percent  Confidence  Limit  •  Observed  B/ents  (Weibull  plotting  positions) 


SWE  in  Inches  ^  SWE  in  Inches 
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Figure  21.  A  frequency  plot  for  the  Old  Enterprise  Mill  snow  course  (OEM). 

General  Frequency  Analytical  Rot  for  OEM 
Return  Period 


—  5  Percent  Confidence  Limit  - 95  Percent  Confidence  Limit 

Observed  Events  (Wei bull  plotting  positions)  □  Low  Outlier 


Figure  22.  A  frequency  plot  for  the  Quaking  Aspen  snow  course  (QKA). 

General  Frequency  Analytical  Rot  for  QKA 
Return  Period 


- 5  Percent  Confidence  Limit  - 95  Percent  Confidence  Limit 

O  Observed  Events  (Wei bull  plotting  positions)  □  Low  Outlier 
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Figure  23.  A  comparison  of  the  annual  probability  estimated  at  each  of  the  SWE  sites  and  the 

initial  SWE  distribution  of  the  PMF. 


Elevation  (feet) 


Next,  we  determined  the  impact  of  the  initial  SWE  distribution  on  the 
snowmelt  calculated  during  the  PMF  analysis.  Figure  24  shows  the  SWE 
that  was  melted  for  each  elevation  band.  One  can  see  that  the  SWE  distri¬ 
bution  can  be  divided  into  four  regions,  starting  at  the  highest  elevations. 

1.  In  the  highest  elevation,  the  snow  did  not  melt  but  actually  increased  dur¬ 
ing  the  PMF  analysis.  This  reflects  the  subfreezing  temperatures  that  ex¬ 
isted  at  this  high  elevation  during  the  event. 

2.  In  the  next  region,  the  SWE  was  greater  than  the  snowmelt  that  occurred 
during  the  PMF.  The  snowmelt  in  this  region  was  heat  transfer  limited. 
The  snowmelt  runoff  during  the  PMF  was  insensitive  to  the  amount  of 
SWE  given  that  the  SWE  was  greater  than  the  snowmelt  amount.  This  re¬ 
gion  covered  the  5000-6000, 6000-7000, 7000-8000,  and  8000- 
9000  ft  elevation  bands. 

3.  In  the  third  region,  the  snowmelt  was  limited  by  the  initial  SWE  amount. 
Runoff  from  this  region  was  sensitive  to  the  amount  of  initial  SWE.  This 
region  covered  the  3000-4000  ft  band  (4  in.  of  initial  SWE)  and  the 
4000-5000  ft  band  (8  in.  of  initial  SWE) 
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4.  The  fourth  region  is  below  the  snow  line  and  contributed  no  snowmelt 
during  the  PMF.  This  region  covered  the  550-1000, 1000-2000,  and 
2000-3000  ft  bands. 

Above  5000  ft,  the  initial  SWE  distribution  was  greater  than  the  snowmelt 
that  occurred  during  the  PMF.  Therefore,  in  this  region,  the  PMF  results 
were  insensitive  to  the  actual  amount  of  SWE  and  the  initial  SWE  distribu¬ 
tion  appropriate.  Below  5000  ft,  the  snowmelt  equaled  the  initial  SWE  dis¬ 
tribution,  indicating  that  the  snowmelt  was  limited  by  the  initial  SWE.  Un¬ 
fortunately  there  is  no  information  on  the  actual  SWE  in  this  region  of  the 
Success  Dam  watershed.  Based  on  the  work  of  others  (Rice  et  al.  2011)  and 
the  SCA  analysis  above,  the  distribution  selected  in  this  region  seems  ap¬ 
propriate. 


Figure  24.  The  calculated  SWE  melted  during  the  PMF  analysis  and  the  initial  SWE  for  each 
elevation  band.  (Note  that  SWE  melted  is  shown  as  positive;  SWE  increases  during  the  PMF 

analysis  are  shown  as  negative.) 
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7  Snow  Melt  Modeling 

The  inclusion  of  SCA  in  the  snowmelt  modeling  will  have  two  major  im¬ 
pacts.  The  first  is  that  the  rate  of  snowmelt  reaching  the  soil  surface  of 
each  elevation  band  will  be  reduced  by  the  SCA  fraction.  This  is  apparent 
when  we  describe  the  rate  of  snowmelt  generated,  R,  as 

R  =  mfA(T.-T,)  (2) 


where 

m/  =  the  melt  factor, 

A  =  the  area  measure  of  the  elevation  band, 

Ta  =  the  air  temperature,  and 

Tm  =  the  melt  temperature  (generally  set  as  32°F). 

When  SCA  is  included,  R  becomes 

R=CSCAmfA(Ta-T.)  (3) 

where  Csca  =  SCA  of  the  elevation  band. 

The  second  impact  is  that  the  total  volume  of  runoff  from  an  elevation 
band  will  be  reduced  by  the  SCA  fraction.  The  impacts  on  the  total  volume 
can  be  significant  as  described  in  the  previous  section.  Snow  models  that 
rely  on  calibration  for  estimating  melt  factors  may  need  to  recalibrate 
when  SCA  is  included. 
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8  Summary 

This  report  reviewed  the  initial  snowpack  distribution  assumed  to  be  in 
place  at  the  start  of  the  PMF  analysis  of  Success  Dam  (USACE  2013).  The 
PMF  snowmelt  analysis  divided  the  watershed  into  elevation  bands  cover¬ 
ing  1000  ft  in  elevation  and  assigned  a  snow-water  equivalent  (SWE) 
depth  to  each  band.  The  initial  SWE  ranged  from  a  maximum  54.4  in.  at 
8000-9000  ft  to  a  minimum  of  4  in.  at  3000-4000  ft.  No  snow  was  as¬ 
sumed  below  3000  ft.  Each  elevation  band  was  implicitly  assumed  to  be 
100%  snow  covered.  The  initial  distribution  of  the  snowpack  resulted  in 
232.7  square  miles  (59.7%)  of  the  watershed  covered  with  snow  and  a  total 
SWE  volume  of  267  x  io3  acre-ft. 

CRREL  analyzed  satellite  imagery  to  determine  the  SCA  of  the  Success 
Dam  watershed  from  2000  to  2014.  This  analysis  was  based  on  TDF  of 
NASA’s  daily  MOD10A1  snow  product.  TDF  results  in  virtually  cloud-free 
daily  snow  maps  displaying  the  fSCA  of  pixels  not  obscured  by  clouds.  The 
pixel  size  is  roughly  500  x  500  m.  The  Success  Dam  watershed  is  covered 
by  4710  pixels.  The  fSCA  was  estimated  for  each  pixel  for  each  day  for 
2000-2014. 

The  CRREL  SCA  analysis  found  that  the  SCA  of  the  Success  Dam  water¬ 
shed  is  highly  dynamic  with  large  day-to-day  variations.  We  considered 
various  options  for  selecting  the  appropriate  SCA  to  use  in  determining  the 
initial  SWE  distribution  and  calculated  the  SCA  statistics  for  each  day  of 
the  winter  season.  We  selected  the  maximum  95th  percentile  SCA  for  each 
elevation  band  was  selected  as  being  sufficiently  conservative  without  be¬ 
ing  extreme.  The  total  watershed  SCA  that  results  from  selecting  the  maxi¬ 
mum  95th  percentile  SCA  for  each  elevation  band  is  101.6  square  miles 
(26.03%  of  the  watershed).  This  is  greater  than  the  total  watershed  SCA 
observed  on  99.5%  of  the  days  with  SCA  greater  than  zero.  The  SWE  vol¬ 
ume  that  results  from  combining  the  SWE  depth  for  each  elevation  band 
used  in  the  PMF  study  with  the  95th  percentile  SCA  for  each  elevation 
band  is  178  x  ios  acre-ft. 

We  reviewed  the  initial  SWE  distribution  used  in  the  PMF  study.  Above 
5000  ft,  the  initial  SWE  distribution  was  greater  than  the  snowmelt  that 
occurred  during  the  PMF.  The  snowmelt  in  this  region  was  therefore  heat 
transfer  limited.  In  this  region,  the  snowmelt  runoff  during  the  PMF  was 
insensitive  to  the  amount  of  SWE  given  that  the  SWE  was  greater  than  the 
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snowmelt  amount.  Below  5000  ft,  the  snowmelt  was  limited  by  the  initial 
SWE.  Runoff  from  this  region  was  sensitive  to  the  amount  of  initial  SWE. 
Unfortunately  there  is  no  information  on  the  actual  SWE  in  this  region  of 
the  Success  Dam  watershed.  Based  on  the  work  of  others  and  the  SCA 
analysis  described  in  this  report,  the  distribution  selected  in  this  region 
seems  appropriate. 

The  100%  snow  cover  area  used  in  the  PMF  study  is  also  a  requirement  of 
the  HEC-HMS  elevation  band  snow  model.  HEC-HMS  is  not  currently  ca¬ 
pable  of  including  the  impact  of  SCA  in  the  elevation-band  snow-model 
option  that  was  used  for  the  PMF  runoff  calculations.  An  alternative  is  to 
use  the  gridded  snow  model  option  instead.  Additionally,  other  ap¬ 
proaches  could  be  used  that  could  conserve  the  volume  of  runoff. 
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